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CAVITY  EIGENMODES  FOR  THE  NIST/NRL  FREE  ELECTRON  LASER 


I.  INTRODUCTION 


The  NIST( formerly  NBS)/NRL  FEL  oscillator  powered  by  a  CV  185  MeV 

1-3 

racetrack  microtron  is  currently  under  construction  .  A  simplified 

model  of  the  resonator  cavity,  consisting  of  the  wiggler  vacuum  chamber 

and  the  two  mirrors,  is  shown  schematically  in  Fig.  1.  Since  stimulated 

emission  takes  place  predominantly  along  the  electron  beam  path,  the 

vector  potential  of  the  radiation  is  expressed  in  terms  of  the  free 

space  eigenmodes  A„(r)  =  A  (r)  exp(-iwt)  e _  +  cc  of  the  paraxial  wave 

k  mp  up 

4 

equation  ,  where  e  is  the  polarization  vector, 


*  ,r) .  VriU)  eik(2  *  e*V2> 

Amp(r)  “  21/2  e 


(1) 


k  =  w/c  is  the  wave  number  and  w  is  the  frequency.  In  Eq.  (1)  the  exponent 
2  2 

k(z  +  (x  +y  )/2R(z)]  contains  the  phase  variation  on  the  wavelength  scale 

2  2 

X  =  2n/k,  with  spherical  wavefronts  of  curvature  1/R(z)  =  z/(z  +  b  ).  The 

slow  phase  variation  is  given  by  C  (z)  =  (2m  +  p  +  1)  tan~*(z/b).  The 

^  22  1/2 

spot  size  of  the  radiation  envelope  is  V (z)  =  w  (1  +  z  /b  )  ,  where  the 

1/2 

distance  z  is  measured  from  the  position  of  the  waist  w  =  (2b/k)  .  The 

amplitude  squared  of  the  mode  drops  by  1/2  over  a  distance  equal  to  the 

Rayleigh  length  b  (also  known  as  confocal  parameter).  Most  of  the 

radiation  is  confined  within  a  cone  parametrized  by  the  dif  a;. ion  angle 

1/2 

=  W/z  a  (X/bn)  .  For  given  wavelength  X,  any  two  the  four 

parameters  R,  W,  b,  z  determine  a  mode  uniquely. 

The  amplitude  profile  ump(r*V)  contains  the  transverse  spatial 

variation,  equivalent  to  a  small  k^,  perpendicular  to  the  z-direction.  In 

?  2  1/2 

cylindrical  coordinates  (r,0,z),  where  tan©  =  x/y,  r  =  (x  +y  )  , 

u  (r,0;V)  takes  the  form 
mp 
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(2) 


u_  (r, 0; V)  =  a 
mp  ’  '  mp 


'cosp©' 

>sinp0. 


lp/2  LfcO  e  2  , 


2r 

=  2  ’ 


where  +  (-)  signifies  cosine  (sine)  poloidal  dependence,  lP(£.)  are  the 
associated  Laguerre  polynomials  and  the  normalizing  factor  amp  = 
ll/(USp0)f.W2]1/2[m!/(m  +  p)!]-172. 

In  the  presence  of  conducting  walls  the  vacuum  expansion  (1)  still 
provides  the  best  representation,  because  (a)  the  chamber  inner  radius  is 
much  larger  than  the  radiation  spot  size  and  (b)  the  transit  time  through 
the  cavity  and  the  length  of  the  light  pulse  itself  are  too  short  tc  allow 
multiple  reflections  on  the  walls  and  set  up  cavity  eigenmodes.  It  can 
also  be  argued  that  the  small  fraction  of  the  radiation  reflected  from  the 
wall,  acting  as  a  perfect  conductor  for  grazing  incidence,  is  lost  out  of 
the  resonator.  The  main  effect  of  the  chamber,  therefore,  comes  from  the 
two  edge  apertures,  where  the  radiation  spot  size  is  maximum.  The 
resonator  is  then  modeled  by  a  sequence  of  four  optical  elements,  i.e.,  two 
apertures  and  two  mirrors. 

The  radiation  profile  is  altered  after  each  encounter  with  an  optical 
element.  A  pure  incident  mode  Amp(r)  will,  in  general,  be  partially 
transformed  into  different  modes.  This  is  caused  by  the  finite  size  of  the 
apertures,  and,  in  addition,  by  spherical  aberration  and  surface 


imperfections  in  case  of  the  mirrors.  Consider  the  incoming  radiation  to  a 
given  optical  element  as  consisting  of  various  modes  (m,p)  of  the  same 
curvature  R^(z).  Both  incident  and  reflected  radiation  are  expanded  into 
eigenmodes,  respectively  as  follows, 


Al(ri>  ■  Vri>  ■ 


A°(r  )  =  V  c°  A  (r  )  .  (3) 

'  o'  "  nq  nqv  o'  v  ' 

n,q 

The  relation  among  the  incident  and  reflected  expansion  coefficients  c1 
.  o 

and  c_  is  written  as 


o  „  1 

c  =  R  c  , 


c°  =  V  Rmp  c1 
nq  ^  nq  mp  ’ 


where  Rmp  are  the  elements  of  the  reflection  matrix  R. 
nq 

The  radiation  profile  at  the  end  of  the  round  trip  inside  the 
resonator  will  relate  to  the  original  profile  through  the  resonator 
transfer  matrix  M, 


M  ~  T2  R2  T2  T1  R1  Tl’ 


where  T^  and  R^  are  the  transmission  matrices  through  the  chamber  apertures 
and  the  reflection  matrices  from  the  mirrors  1  and  2  respectively.  The 
cavity  eigenmodes  Cj  with  eigenvalues  \k  are  given  by 


M  C.  =  v.  C. 
J  J  J 
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Since  M  is  generally  nondiagonal,  the  eigenmodes  are  mixtures  of  vacuum 
modes  ( 1 ) . 


The  electron  beam  is  an  active  medium  that  changes  the  radiation 
profile  during  amplification  in  each  passage.  If  G  is  the  amplification 
matrix,  then  a  steady  state  exists  finally  if  the  matrix  equation 

G  M  C  =  g  C  (7) 

has  solutions  with  |gg  |  =  1.  A  steady  state  need  not  be  an  eigenmode  of 
the  empty  resonator;  this  could  happen  only  in  the  case  of  equal 
amplification  g  for  each  eigenmode,  i.e.,  G  =  g  I  where  I  is  the  identity 

SL  '"1 

matrix.  In  an  FEL,  a  different  gain  is  associated  with  each  free  space 
eigenmode.  However,  in  cases  when  the  off-diagonal  elements  of  both  G  and 
M  become  vanishingly  small,  both  the  resonator  modes  and  the  final  steady 
states  approach  the  pure  vacuum  modes  (1). 

In  this  paper  we  first  study  the  vacuum  performance  of  the  optical 
cavity.  The  detailed  cavity  mode  structure  in  terms  of  vacuum  modes  and  the 
associated  eigenvalues  are  obtained  by  numerical  diagonalization  of  the 
cavity  matrix  M.  The  fractional  power  loss  per  cavity  mode  per  trip  is 
found  from  the  magnitude  of  the  eigenvalue 

<8, 

The  profiles  for  the  cavity  modes  are  also  obtained  utilizing  the  expansion 
coefficients  of  into  the  vacuum  modes.  We  also  examine  the  eigenmode 
structure  of  the  combined  gain-transfer  matrix  G  M  in  cases  of  small  gain. 


A 


II.  TRANSFER  MATRIX  FOR  A  SINGLE  OPTICAL  ELEMENT 
The  reflection  of  Gaussian  light  beams  from  mirrors  was  studied  in 
some  detail  in  Ref.  5  for  arbitrary  angle  of  radiation  incidence.  In  the 
limit  of  normal  incidence  considered  here,  the  reflection  matrix  elements 
are  giver  by  the  surface  integrals 


lmp 

nq 


d0  dr  r 
s  s  s 


u  (r,0)u  (r,0) 

mpv  s’  s'  nqv  s’  s 


2  ,1/2 


,1/2 


1  + 


z  (r  ) 
s  v  s' 


1  + 


z  (r  ) 
s  v  s 


b.‘ 

i 


iA(rs'0s) 


i 

x  e 


(9) 


The  mirror  surface  S  is  spherical,  expressed  in  the  coordinate  system 
(rs,0s,zs)  with  origin  located  at  the  mirror  center,  by 


(10) 


where  R  is  the  mirror  radius  of  curvature.  Equation  (10)  is  used  to 
m 

express  z  on  the  surface  S  in  terms  of  r  .  The  mirror  boundary  is  given 
s  s 

by 


(11) 


where  p  is  the  radius  of  the  mirror  cross-section. 

The  phase  factor  A(r  , ©  ),  related  to  the  optical  path  along  the 

s  s 

various  rays  connecting  the  incoming  wavefront  with  its  mirror  image 
(reflected),  must  be  approximately  constant.  Therefore,  the  curvature  of 


5 


the  outgoing  wavefront  is  related  to  the  incoming  and  the  mirror  curvatures 
through 


(12) 


The  Rayleigh  length  bQ  and  the  waist  location  1  of  the  outgoing  modes  are 
yet  to  be  determined. 

It  has  been  argued^’^  that  the  amount  of  radiation  scattered  into 
other  than  the  incoming  modes,  as  well  as  that  escaping  behind  the  mirror, 
depends  on  three  factors: 

(a)  Finite  mirror  size  effects,  of  the  order  of 


where  p  is  the  mirror  cross-section  radius  and  the  incoming  radiation 
spot  size. 

(b)  Spherical  aberration  effects,  of  the  order  of 


,  w.  , 

?i  (  71  )  ’ 


coming  from  the  phase  term  a(r  ,0  )  inside  (9).  Spherical  aberration 

s  s 

exists  even  when  the  mirror  curvature  matches  the  radiation  curvature.  It 
is  caused  by  the  fact  that  rays  ending  on  a  given  spherical  wavefront  are 
not  exactly  perpendicular  to  it,  since  they  originated  from  a  finite  size 
waist  and  not  from  a  point  at  the  center  of  curvature. 

(c)  Surface  imperfections,  for  example,  when  the  reflecting  surface 


is  not  perfectly  spherical. 

Spherical  aberration  effects  are  usually  less  important.  They  will 
be  addressed  in  future  work,  together  with  the  potentially  more  important 


consequences  caused  by  mirror  deformations  (buckling)  due  to  heating. 


6 


Ignoring  (b)  and  (c)  amounts  to  setting 
substituting  expression  (2)  for  the  eigenmodes  u 


A(r 

mp’ 


g,0  )  =  0  in  (9).  After 
(9)  becomes, 


d(u’°) 


=  5 


-.mp 
'pq  ~nq 


d£  (a202 


E 

rl 


9/^2 , 


( a  £)  1/(0 

n  m 


(15) 


0 


where  cj^  =  a  [m!n!/(m+p) ! (n+q) \ J*^2  and  £  =  2r2/V^2.  Since  the  surface  S 
has  rotational  symmetry  about  the  2-axis,  it  couples  modes  with  the  same 
poloidal  9  dependence,  p  =  q.  The  radial  integration  is  carried  out  in 
Appendix  A. 

In  general,  the  matrix  R  involves  two  independent  parameters,  the 

ratio  of  the  mirror  radius  p  to  the  radiation  spot  size  squared, 

2 

y  =  (p/Wj)  ,  and  the  ratio  of  the  incoming  to  outgoing  spot  sizes 
a  =  (Vj/Wq).  Only  the  curvature  of  the  reflected  mode  is  set  by  the 
mirror,  while  the  outgoing  spot  size  is  still  a  free  parameter.  This  can 
be  exploited  by  choosing  the  value  Wq  that  maximizes  the  coefficient  for 
the  fundamental  mode  in  the  reflected  radiation,  i.e., 


3R 


00 

00 


=  0  . 


3a 


(16) 


Once  Wq  is  selected,  the  exact  location  and  size  of  the  waist(s)  for  the 
reflected  modes  is  determined  by  solving  the  system  of  equations 


V 


o 


w 


o 


1/2 

9 


1 

R 


l2+b2 
o  o 


(17) 


The  transmission  matrix  T  through  an  aperture  is  given  by 


T  = 


(18) 
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where  R  is  the  reflection  matrix  for  a  plane  mirror  (R  =  ®)  of  the  same 

m 

cross  section  p,  while  the  curvature  transformation  is 


(19) 


Aberration  and  surface  imperfections  do  not  affect  transmission  through 
apertures. 


8 


III.  CAVITY  EIGENMODES 


Ve  seek  a  class  of  cavity  igenmodes  with  the  waist  located  in  the 

middle  of  the  vacuum  chamber.  Since  the  mode  spot  must  remain  unchanged 

V.  =  V  during  each  transmission  or  reflection  we  elect  a  =  1  inside  the 
10 

transfer  matrix  (5)  for  every  optical  element.  The  curvature  R  is  not 
changed  during  transmission  through  an  aperture.  The  mirror  curvature, 
however,  must  match  the  incoming  and  outgoing  radiation  curvature,  R^  =  Rq 
=  Rm  in  (12).  Therefore,  for  eigenmodes  to  exist,  the  equations 


R1  =  2  2  ' 
1  +  bz 


R2  =  .2  .  2  ’ 
+  b 


must  admit  a  positive  solution  for  b.  This  is  possible  when 


L  -  L1  +  L2  * 


Equation  (21)  is  the  optical  stability  condition  for  the  cavity.  Even 

with  (21)  satisfied  the  cavity  modes  will  still  decay  slowly  in  time 

because  of  the  finite  size  of  the  mirrors.  According  to  (20),  the 

Rayleigh  length  b  for  various  wavelengths  remains  fixed  for  a  given  mirror 

1/2 

configuration,  while  the  waist  w  varies  as  w  =  (bA/it)  .  The  cross 

section  of  the  electron  beam  is  adjustable  to  give  a  good  filling  factor 

for  the  particular  wavelength.  The  cross  section  of  the  chamber  is  oval 

1/2 

and  is  approximated  with  a  circular  one  of  effective  radius  p  =  (p^^)  > 

where  p^  and  P2  are  the  major  and  minor  aperture  radii. 

Because  of  the  imposed  axisymmetry,  we  look  for  cavity  eigenmodes 


involving  combinations  of  axisymmetric  vacuum  modes  p  =  q  =  0.  The 


coefficients  c.  '  of  the  modes  u  n  inside  the  j-th  cavity  mode  C. 

J  nU  ^  j 


C. 

J 


(22) 
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and  the  corresponding  eigenvalues  are  found  by  numerical 
diagonalization  of  the  transfer  matrix  H,  utilizing  expressions  (A3)  with 
a=  1  for  each  individual  optical  element  inside  (5).  Eventually,  the 
cavity  matrix  M  depends  on  the  spot  to  aperture  size  ratios  for  the  four 
optical  elements, 

MIu^X),  p2(X),  y3(X),  u4(X)]  ,  (23) 


where  in  turn  y^  depends  on  the  wavelength  through 


(24) 


Two  general  conclusions  are  made.  First,  it  is  found  that,  for 

radiation  spot  sizes  smaller  than  about  one  third  the  aperture  sizes,  the 

dominant  contribution  in  each  cavity  eigenmode  comes  from  a  single 

coefficient  c/n^.  Cavity  eigenmodes,  in  this  case,  approach  pure  vacuum 

modes,  as  expected  from  the  smallness  of  the  off-diagonal  matrix  elements 

in  M.  Second,  if  the  eigenvalues  are  arranged  according  to  magnitude, 

|\>0|  >  |vj|  >  ...  |Vj  |  >  lvj  +  il  the  largest  eigenvalue  corresponds  to 

the  eigenmode  closest  to  the  fundamental  vacuum  mode  uqq>  i*e*»  the 

eigenmode  with  |cq^  |  «  1.  The  next  largest  eigenvalue  corresponds  to 

the  eigenmode  closest  to  the  first  radial  vacuum  mode,  i.e.,  with  |c^^  | 

=  1,  and  so  on.  This  is  expected  as  the  rms  spot  size  for  the  n-th  radial 

1/2 

mode  unQ(r;V)  increases  with  n  as  (n+1)  W. 

The  NIST/NRL  oscillator  has  been  designed  for  eigenmodes  of  Rayleigh 
length  equal  to  half  the  vacuum  chamber  length.  Two  different 
arrangements,  one  with  full  ’•  -r-;ler  and  one  with  half  length  wiggler  will 
be  used  for  the  wavelengt  jgimes  of  0.2ym  to  2ym  and  2ym  to  lOym 
respectively.  The  design  para"  '  rs  are  b  =  IW2  =  107.5cm,  Lj  =  521cm, 


10 


L2  =  386cm,  R^  =  543cm,  Rj  =  417cm,  ~  2.54cm  for  the  half  viggler 
and  b  =  Ly/2  =  198cm,  =  431cm,  =  477cm,  R^  =  521cm,  R2  =  559cm  for 
the  full  viggler. 

2 

The  fractional  power  loss  hj  =  1  -  |\k  |  for  the  first  5  cavity 
eigenmodes  of  the  half  viggler  arrangement  is  plotted  in  Fig,  2  as  a 
function  of  the  wavelength  X.  Through  the  planned  regime  of  operation  the 
loss  for  the  fundamental  cavity  mode  never  exceeds  IX.  The  loss  factor 
for  the  next  two  modes  is  also  very  small,  so  that  the  mode  selection  is 
going  to  be  determined  by  the  differences  in  the  radiation  gain  for  each 
mode.  In  Fig.  3(a)  we  plot  the  expansion  coefficients  c  of  the 
fundamental  cavity  mode  Cq  into  the  vacuum  modes  unQ.  It  shows  that  the 
cavity  mode  profile  is  very  close  to  a  pure  Uqq  vacuum  mode,  with  other 
modes  contributing  less  than  1%.  The  expansion  coefficients  for  the 
second  and  third  cavity  modes  are  shown  in  Figs.  3(b)  and  3(c) 
respectively.  If  the  complete  expansion  coefficients  are  written  as 

cj  =  lcj  I  exP  <ixjn)  , 


then  the  complex  amplitude  for  the  j-th  eigenmode  Cj  is  given  by 


i*.(r) 

Aj(r)  =  lAj(r> I  e  3 


where 


1/2 

I A j ( r )  |  =  [  A r^2  +  AijZ  ]  ,  $j(r)  =  tan-1  [  Ar^  /  Aij  ]  , 


Arj  -  2  »n0<r)  lcj<n)|  cos  Xjn 


<n)|  ™  v  ,  Aij  =  E  un0(r)  lcj<n>l  sin  xjn 
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The  resulting  amplitude  profiles  for  the  first  three  modes  at  z  =  0  are 

shown  in  Fig.  4  for  X  =  2pm  and  X  =  10pm.  In  Fig.  5(a)  we  show  the  power 

fraction  t|-  for  the  first  four  cavity  modes  of  Rayleigh  length  b  =  L  /6. 

3  v 

The  vacuum  expansion  coefficients  Cq^1^  f°r  the  fundamental  cavity  mode 
are  shown  in  Fig.  5(b)  against  the  wavelength  X.  The  transverse  amplitude 
profiles  for  the  first  three  modes  at  X  =  2.2pm  and  X  =  10pm  respectively 
in  Figs.  5(c)  and  5(d),  show  considerable  departure  from  vacuum  modes. 

The  cavity  eigenmodes  for  the  full  wiggler  arrangement  are  extremely 
close  to  vacuum  modes  and  the  fractional  power  losses  for  the  first  five 

_3 

of  them  are  below  10  over  the  frequency  regime  from  0.5pm  to  2pm.  This 
is  caused  by  the  combination  of  a  longer  Rayleigh  length  with  a  spot  size 
that  gets  smaller  with  shorter  wavelength. 

In  case  of  small  gain  we  can  include  the  effect  of  the  electron  beam 
on  the  cavity  eigenmodes  by  introducing  the  amplitude  gain  matrix  G.  In 
the  linear  regime  the  cross-coupling  among  various  transverse  modes  is 
unimportant  and  G  is  diagonal,  given  by 


G  =  g(X)  f  6  , 

mn  6  v  '  n  mn 


(25a) 


where 


g(X)  =  0.5  Fj2 


(25b) 


is  the  amplitude  gain  for  the  fundamental  vacuum  mode.  In  Eqs.  (25)  N  is 
the  number  of  wiggler  periods,  r  is  the  initial  relativistic  factor,  = 

O  K 

2  3 

itw  is  the  radiation  cross  section  at  the  waist,  1^  =  17x10  A,  I  is  the 

2 

current  in  amperes,  K  =  |e |B  X  /2nmc  is  the  wiggler  parameter,  B  is  the 

V  V  V 

rms  magnetic  field  of  the  wiggler,  Xy  is  the  wiggler  wavelength,  related 
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9  2 

to  the  radiation  wavelength  X  by  Xy  =  2y‘'X/(l+K  ),  JQ(b)-Jj(b)  with  b 
2  2 

=  K  /2(1+K  )  and  f  is  the  normalized  (to  the  fundamental)  filling  factor 


dr  r  j(r)  u 


nO 


<r)  /  I 


dr  r  j (r) 


u00(r) 


with  the  parabolic  current  profile  given  by  j(r)  =  j 
pjj»  j(r)  =  0  for  r  >  p^,  being  the  beam  radius, 

for  the  j-th  beam-cavity  eigenmode  is  given  by 


(26) 

(1  -  r^/p^)  for  r  < 
The  round  trip  gain 


gj  =  i  i  +  g<x)  i 


(27) 


where  is  the  eigenvalue  for 
transfer  matrix 


the  eigenmode 


cj 


of  the  combined  gain- 


M  =  G  M  .  (28) 

g 

The  fractional  power  loss  hj  =  1  -  |  for  the  first  five  eigenmodes  is 

shown  in  Fig.  6  as  a  function  of  the  wavelength,  while  the  expansion 
coefficients  of  the  fundamental  eigenmode  in  terms  of  vacuum  modes  are 
shown  in  Figs.  7(a)-7(c).  The  losses  for  the  higher  eigenmodes  j  >  0  are 
now  considerably  higher  than  the  fundamental  mode  (compare  Figs.  2  and  6), 
thusly,  mode  selection  among  transverse  modes  occurs  through  amplification 
of  the  radiation,  because  of  the  differences  in  the  filling  factor  fR. 
The  transverse  amplitude  profiles  at  X  =  2.2  ym  and  X  =  10  ym  appear  in 

Figs.  8(a)  and  8(b). 
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APPENDIX  A.  Computation  of  the  Reflection  Matrix  Element 


The  associated  Laguerre  polynomials  are  given  by 


m 


Lfa)  =  E  (-D 


(m+p)! 


*k  • 


(Al) 


k=0  k!  (m-k) !  (p+k) ! 


Substituting  (Al)  inside  (15)  and  integrating  by  factors,  using 

m  m! 


w 


n  0-5 

e  =  -  e 


►m-k 


(m-k) ! 


]■ 


(A2) 


one  obtains 


=  6  ^  |m!n!(m+p)!(n+q)! j 


1/2 


+1 


m  n  ,  2  Nk+1 

T  E  E  *2k  — 

1  n  '  n  V +1J 


p+a^Jp+i  k=0  1=0 


(~l)k+1(k+l+P) ! 


k!l! (m-k) ! (n-1) ! (p+k) ! (p+1) ! 


1  -  e-(“2+1)y  kE+P(2ym)k+1+p 
i=0 


,  2  iiP-i 


(k+l+p-l) ! ) 


(A3) 
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Fig.  2  —  Fractional  per  round  trip  power  loss  versus  wavelength  for  the  first  five  resonator  eigenmodes  in 
vacuum.  The  Rayleigh  length  b  is  half  the  wiggler  length.  The  cavity  parameters  for  the  half  wiggler  arrange¬ 
ment  are  L|  =  521cm,  L2  =  386cm,  Rx  =  543cm,  R2  =  417cm  and  p,  =  p2  =  2.54cm. 
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^(A*m) 

Fig.  3a  —  Expansion  coefficients  into  vacuum  modes  for  the  fundamental  cavity  eigenmode  versus  wavelength. 
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Fig.  3b  —  Expansion  coefficients  into  vacuum  modes  for  the  first  cavity  eigenmode  versus  wavelength. 
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Fig.  5b  —  Plots  of  expansion  coefficients  of  the  fundamental  cavity  mode  in  vacuum  modes  versus  wavelength 
where  the  fundamental  cavity  mode  structure  has  Rayleigh  length  of  b  =  Lw/6  and  corresponding  mirror  curva¬ 
tures  of  /?i  =  523cm  and  R2  =  389cm. 


2- 


2.0 


Fig.  5c  Plots  of  transverse  profiles  for  the  first  three  cavity  modes  at  wavelength  2.2  /xm  where  the  funda¬ 
mental  cavity  mode  structure  has  Rayleigh  length  of  b  =  Lw/ 6  and  corresponding  mirror  curvatures  of 
R  i  =  523cm  and  R2  =  389cm. 
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27 


Fig.  7a  —  Structure  of  the  combined  beam-cavity  eigenmodes.  Same  notation  as  in  Fig.  3. 
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^(/i  m) 


Fig.  7b  —  Structure  of  the  combined  beam-cavity  eigenmodes.  Same  notation  as  in  Fig.  3. 
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Fig.  8a  —  Profiles  of  the  first  three  beam-cavity  eigemnodes.  Same  notation  as  in  Fig.  4. 
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